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Spectrum and decay matrix elements of B and D-mesons in lattice 
NRQCD 

Joachim Hein*^* 

'^Newman Laboratory of Nuclear Studies, Cornell University, Ithaca, NY 14853, USA 

We discuss recent results on the excitation spectra of B and D-mesons obtained in the framework of non- 
relativistic lattice QCD in the quenched approximation. The results allow for the determination of the MS-mass 
°^ "^6 MS^''^b Ms) ~ 4.34(7) GeV in 0{as) in the perturbative matching. The determination of the decay constants 
/sg and /oa is discussed in detail. Results for the matrix elements of semi-leptonic B to D decays are shown. 



1. MOTIVATION 

In the standard model of elementary particle 
physics, CP-violation is caused by a single phase 
in the CKM-matrix. At the present time it is ex- 
perimentally observed only within the Kaon sys- 
tem. Establishing CP-violation also in the B- 
meson system is one of the major goals for the 
B-factory experiments, such as the newly running 
BaBar, Belle and CLEO-III experiments as well 
as the future hadron collider experiment LHCb. 
Failure of all CP-violating processes to be de- 
scribed by the single CKM-phase would provide 
direct evidence for new physics beyond the stan- 
dard model. 

The extraction of the elements of the CKM- 
matrix from the above experiments is complicated 
due to the hadronic nature of the initial and fi- 
nal states. Here accurate knowledge of hadronic 
quantities, such as QCD form factors of the B- 
meson, is needed. Due to confinement these quan- 
tities are genuine non-perturbative. Lattice gauge 
theory provides a means to determine such prop- 
erties of hadronic states based on first principles 
in a model independent way. The approximations 
made in present day calculations are expected to 
be relaxed in future calculations. 

2. LATTICE NRQCD 

For the results presented here we use non- 
relativistic QCD (NRQCD) ||i|,|| to formulate 



heavy b and c-quarks on a lattice with a lat- 
tice spacing a which is not negligible against its 
Compton wave length. For heavy-light mesons, 
this results in a systematic expansion of the 
Hamiltonian in powers of the inverse heavy quark 



H = Ha + 6H + SH, 
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The leading kinetic term is given as 
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where D denotes a covariant lattice derivative. 
We also include the following relativistic correc- 
tions 



5H ■ 



9 



2toc 

g 

'8m| 



8m%, 



a-(DxE-ExT>)- 



8m%, 



(3) 



To reduce the dependence on the lattice spacing, 
the following discretisation corrections are also in- 
cluded 
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The first term on the right-hand side corrects in 
spatial, the second one in temporal direction. The 
n denotes the stability parameter used in the evo- 
lution equation to generate the heavy quark prop- 
agators. The inclusion of the improvement terms 



is important due to the non-renormalisability of 
NRQCD. The action has to be improved at finite 
values of the lattice spacing such that the results 
become independent of a within the achieved ac- 
curacy. We did not include radiative corrections 
to the prefactors of the individual terms of this ex- 
pansion. However, the dominant tadpole contri- 
butions to this coefficients are removed by mean 
field theory g]. 

The light quarks are simulated with the clover 
action El and the standard plaquette action is 
used for the gluon background. All the results 
reported here are calculated in the quenched 
approximation. This means the neglection of 
vaccum polarisation effects due to the light sea 
quarks. For some of the results also partly un- 
quenched results exist. These calculations include 
two flavours of light sea quarks. This will be 
pointed out in the individual subsections. 

In the calculations the value of the lattice spac- 
ing is determined from the mass of the p-meson. 
This is senseful, since the here discussed B and D- 
mesons have a typical momentum scale of Aqcd 
as does the p. In the quenched theory the strong 
coupling constant runs differently from nature. 
This procedure minimises the quenching error 
compared to using a quantity of a different scale. 

3. MESON SPECTROSCOPY 

Due to heavy quark symmetry the excitation 
spectra of the B and D-mesons are highly re- 
lated. The lowest lying states are experimentally 
well understood [p[, however the situation with 
respect to radially and orbitally excited states is 
not as satisfactory. In this situation we can learn 
about the accuracy of the approximations done in 
the present calculation, as well as predict masses 
for so far undiscovered states. 

3.1. Spectrum of the B 

The spectrum of the i3-meson has been in- 
vestigated in large detail in ^,0]. In summary 
these publications provide results for three differ- 
ent values of the lattice spacing a. Within the 
achieved accuracy, the individual results for the 
different splittings agree with each other. This is 
particularly important because of the above dis- 
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Figure 1 . Excitation spectrum of the 5-meson. 



cussed non-renormalisability of NRQCD. A sub- 
set of the investigated states are also reported by 
the JLQCD Collaboration B using similar meth- 
ods. These results are in excellent agreement. 

The results of MM are summarised in figure n^. 
The lattice results are displayed by the octagons, 
whereas experimental results are given by hor- 
izontal lines. We compare the lattice B** and 
B** results to the experimental i3}(5732) and 
B*j{5850) resonances. The dashed line gives a 
preliminary DELPHI result. The overall agree- 
ment between lattice and experiment is very good 
and the spin-independent spectrum is nicely re- 
produced. This includes the radial and orbital 
excitation energies as well as the strange non- 
strange 5- wave splittings. However there are 
problems with the spin-dependent hyperfine split- 
ting, which turns out to be significantly too small. 
This might be due to the neglected radiative cor- 
rections in eqn. (g), especially in front of the crB 
term. Preliminary results P| indicate this might 
have an effect of the order of 10% on the hyperfine 
splitting. Another cause might be the quenched 



approximation. Similar problems are also ob- 
served in quenched light spectroscopy. The hy- 
perfine splitting turns out increasingly too small 
with increasing quark mass |Q . So far results for 
the B-meson hyperfine splitting with the inclu- 
sion of two flavours of light sea quarks do not show 
any improvement for fixed lattice spacing ||ll|,|l2[ . 
However these calculations use quite large values 
for sea quark mass and improvement might only 
be seen once more realistic values are used. 

For the first time, this calculation gives a result 
for the radial excited P-states, which is denoted 
a.s B**{2P) in the figure. 
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3.2. Spectrum of the D 

NRQCD calculations of the i?-meson within 
the here-presented framework are only reliable for 
larger values of the lattice spacing. The conver- 
gence of NRQCD for £)-mesons has been inves- 
tigated in |1^]. A detailed analysis of the pre- 
sented material shows reasonable convergence of 
the NRQCD expansion for charmed heavy light 
states 0. 

The most complete results to date on the D- 
meson excitation spectrum have been published 
in 0. These results, using a lattice spacing of 
0.177 fm, are presented in figure 0. Again the 
symbols give the lattice results and the horizon- 
tal lines the experimental outcome. The dashed 
line gives a result of the DELPHI collaboration 
and the shaded region shows a preliminary re- 
sult of the CLEO collaboration for a wide Di- 
resonance. Again the spin-independent spectrum 
agrees very well with the experimental outcome. 
Also the experimentally observed increase of the 
strange non-strange splitting from the B to the 
I^-system is well reproduced. Again the hyper- 
fine splitting turns out to be too small however 
in this case by not as much as for the B. This 
is consistent with the observed mass-dependence 
of the hyperfine splitting for light quarks in the 
quenched approximation [[lO| . 

Some of the splittings in figure g have also been 



computed with heavy clover quarks 1 14 15f| and no 



significant differences have been found between 
the heavy clover and the NRQCD result. This 
is particularly remarkable for the D* — Dg hy- 
perfine splitting, which strongly depends on the 
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Figure 2. Excitation spectrum of the D-meson 
for a-^ = 1.116 GcV. 



details of the applied heavy quark action. For 
charm quarks, NRQCD and heavy clover quarks 
behave differently. The details of this comparison 
are discussed in 0. 

4. MS-MASS OF THE 6-QUARK 

The mass of the 6-quark is a fundamental pa- 
rameter of the standard model. Based on the 
results of nM it is possible to determine mi, in 
NNNLO perturbation theory. Here preliminary 
results of a forthcoming publication are presented 

4.1. Calculation of the 6-mass 

The calculation is done in two steps. First the 
pole mass is calculated, which is converted into 
the MS-mass in a second step. The calculation of 
the pole mass turns out to be most precise using 
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The mass m-g of the spin-average of the B and 
the B* is taken from experiment [p[ and the sim- 
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Figure 3. Dependence of rn^ Ms ("^ 
lattice spacing a. 
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ulation energy i^sim f^rom lattice simulation. Here 
the use of the spin-averages eliminates the prob- 
lems encountered with the hyperfine splitting in 
the spectrum calculation. The self-energy Eq of 
the 6-quark is known from perturbation theory. 
The latter is available most precise in the limit of 
infinite quark mass. It has been analytically cal- 
culated to 0[a1) Q and numerically to 0{a^g) 
for the quenched cas e |18| , |19| . The here- presented 
analysis is based on [|18[. The simulation energies 
of 1^,0 have been extrapolated to the static limit. 
To convert the pole mass to the MS scheme, we 
use the renormahsation constant to O(a^) po| . 
The result for three different values of the lattice 
spacing is shown in figure 0. The different re- 
sults agree within their statistical accuracy and 
the most precise of them is quoted as the final 
result 



(m,j^)= 4.34(7) GeV. 



"fa.MS^ b,MS 



(6) 



The uncertainty is still dominated by the un- 
certainties of Eq. It further includes the uncer- 
tainties arising from the conversion to the MS- 
scheme, statistics of i?sim, the lattice spacing 



and the discretisation. Corrections due to the 
finite mass of the &-quark are expected to be of 
C(Aqcd/'TI(,) of i?sim, which is about 0.05 GeV. 
The effect of the quenched approximation can 
be estimated at 0{a1) [|7|, when using -Esim with 
two flavours of light quarks from |1^. In com- 
parison to the quenched case in the same or- 
der of perturbation theory, the sea quarks reduce 
"'•b MS ("^b Ms) ^y ^■^'^ GeV. Eqn. (||) is compati- 
ble to the'result in 0{al) of (l|]. 

5. LEPTONIC DECAY 

The pseudo-scalar decay constant /m deter- 
mines the decay of a mesonic state M into a pair 
of leptons 



Pp/m = {Q\A^\M(p)) . 



(7) 



For the _B+-meson the leptonic width is highly 
CKM suppressed and is not expected to be mea- 
surable in the near future. Therefore the decay 
constant cannot be measured directly, however it 
is an important input parameter in other analy- 
sis, for example B — S-mixing ||2l|]. The decay 
constant of the B has to be determined from the- 
ory. The Ds has a much weaker CKM suppression 
and its decay constant has been measured exper- 
imentally, which provides us with the possibility 
to check the theoretical calculation. 

When using NRQCD for the heavy quarks, the 
matrix element of the axial current of QCD has to 
be expanded into matrix elements of the effective 
theory 
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The coefficients Ci on the right-hand side have 
been determined in 1-loop lattice perturbation 
theory |2|]. 

5.1. Investigations of fs^ 

The dependence of the result of /b^ on the lat- 
tice spacing is one of the key issues of [P3i . The 



results have been combined with those of |24| and 
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Figure 4. Lattice spacing dependence of fs. 



Figure 5. Investigation for momentum depen- 
dent discretisation effects on Jb^ ■ 



are reproduced in figure 0. A similar investiga- 
tion has been performed in Iq] . These results are 
included in the figure as well. The individual 
results agree well with each other and one can 
conclude that the results are within the claimed 
accuracy indeed independent of the value of the 
lattice spacing a. For the final number we quote 



Jb^ = 187(16) MeV . 



(9) 



A detailed breakdown of the error can be found in 
23 1 . The result agrees well with the outcome of 
recent studies using other ways to formulate the 
heavy quark on the lattice, see pq,Eq| for reviews. 
Calculations including two flavours of sea quarks 
indicate an increase of approximately 30 MeV for 

Reference [ p3| investigates power law terms and 
the Aqcd/'tiq corrections of the decay constant 
in large detail. Due to the different dimensional- 
ity of the operators on the right hand side of eqn. 
(H) the coefficient Cq develops a l/arriQ diver- 
gence, which cancels a similar divergence in the 
matrix element (^7570(7 • D)*?)- The perturba- 
tive expansion of Co allows for the explicit investi- 
gation of this cancellation in O (as I arriq ) . The re- 



maining part of the ©(I/ttt-q) matrix elements in 
eqn. (0), which contains the 0{p?^laniQ) part of 
the divergence as well as physical 0(Aqcd/"t-q) 
terms, amounts to a few percent of the final re- 
sult. 

The decay constant Jb^ is a momentum inde- 
pendent form factor. It allows the study of the 
effect the lattice discretisation has on the form 
factor of moving mesons. This is an important 
consistency check with respect to semi-leptonic 
decays. In this case the form factors are mo- 
mentum dependent, and discretisation effects and 
physical effects are hard to disentangle. In fig- 
ure H we plot the ratio 
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This has been investigated for two different values 
of a |2^. In the absence of discretisation effects 
this becomes y/ E {p) / E {0) , which is included as 
a full line. The results are consitent with the 
continuum expectation up to momenta of about 
1.2 GeV. For the largest momenta the deviations 
are less than 8%. 



5.2. Decay constant of the Dg 

The decay constant of the Hghter Dg was suc- 
cessfully determined in lattice NRQCD ||2^ as 
well. As for the spectrum, this was done on coarse 
lattices. The result includes all the corrections up 
to ©(Aqcd/'tiq). After cancellation of the power 
divergence in 0{as/amc) the residual 0{l/mc) 
matrix elements contribute only on the 10% level 
to the final number, indicating good convergence 
of the NRQCD expansion in the charm region. 
The final outcome is 



/d, == 223(54) MeV . 



(11) 



This result agrees well with the experimental re- 
suh of fo, = 280(19) (28) (34) MeV H as well 
as other lattice calculations in the quenched ap- 
proximation pa,G6|. 
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6. SEMI-LEPTONIC DECAYS B -^ Dlv 

Semi-leptonic decays of a i?-meson into 
-D, D* , D** , D' , . . . provide the best way to mea- 
sure I Vet I . A precise measurement of this is re- 
quired to relate CP-violation as measured from 
the K to the one measured from the B. In case of 
pseudo-scalar decay products such as D and D' , 
one needs to determine the matrix element 



Form factor /i^ in the case of elastic 



{B\Vf,\D) = y/mBmD[h+{uj){vB+VD)^ 
+h^{ui){vB-VD)fj.] ■ 



(12) 



As shown, this can be parametrised by two form 
factors /i+ (io) and h~{uj) with w = vb'Vd- In this 
section we summarise the results of |g§] , which are 
still preliminary. 

6.1. Elastic scattering 

The elastic scattering of a B-meson from a vec- 
tor current provides the simplest approximation 
to the matrix element in eqn. dl2). In the static 
limit, ViiB — > oo this becomes the Isgur-Wise func- 
tion S,{ijj). The w-dependence of the elastic form 
factor has been studied for two different values of 
the mass of the heavy quark. The result is re- 
produced in figure o. The values for amq — 4 
corresponds approximately to a Bs-nieson. The 
figure gives 



Figure 6 

scattering as a function of the recoil parameter to 



2 
rstranffc 



= 1.5(3)(4) 



(13) 



for the slope of the Isgur-Wise function in case of 
a strange spectator quark. The results have been 
checked for their dependence on the momentum 
of the external states. 

6.2. Radial excited states 

The use of two different interpolating fields 
with the same quantum numbers makes it pos- 
sible to observe a signal for the matrix element 
{B s\VQ{q)\B'^) involving a radial excited state. 
The q dependence of the matrix element is shown 
in figure |7[ For both external states the same 
heavy quark mass is used. The external states 
are orthogonal at zero recoil and the matrix ele- 
ment vanishes as expected. Once these states get 
boosted with respect to each other, the matrix 
element becomes finite. This result is qualita- 
tive and demonstrates the feasibility to determine 
such matrix elements in lattice gauge theory. 

6.3. Non-massdegenerate transitions 

The case of unequal heavy quark mass for the 
in and outgoing state has been studied at zero re- 
coil. In this case the vector current is not a con- 
served current anymore and gets renormalised. 
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Figure 7. Form factor for transitions between 
pseudoscalar ground and radially excited state. 



Figure 8. Form factor /i^ at zero recoil. 



The renormalisation constants have been pertur- 
batively calculated to 0{as) l^]. Figure || gives 
the result with and without this renormalisation 
for three different pairs of heavy quark masses. 
The used heavy quark masses are indicated above 
the individual results in units of the lattice spac- 
ing. Again the value of amq = 4 corresponds to 
the i?s-meson. The arrows at the right hand side 
of the the figure give the upper and lower error 
bound of the lattice result from |30|. This has 
been extrapolated to the physical B ^ D transi- 
tion and corresponds to (l/amb—l/amc)'^ « 0.56. 

7. DISCUSSION 

This talk summarises our recent results on the 
physics of heavy light mesonic states. The results 
include the B and ZJ-meson spectrum as well as 
the decay constants of the Bs and the Dg- The 
results for the B and the Bg have been obtained 
for three different values of the lattice spacing and 
the results are found to be independent of the lat- 
tice spacing within error bars. In the spectrum we 
observe good agreement to experiment for spin 



independent splittings. However the spin depen- 
dent hyperfine splitting comes out too small. This 
has to be resolved in future calculations. 

The spectrum calculations also allow for the 
precise determination of the 6-quark mass in the 
MS-scheme. The preliminary result is m^^ jj^ = 
4.34(7) GeV at its own scale. This result is not 
affected by the too small hyperfine splitting. 

The results for the D-spectrum and the Dg 
decay constant agree well with experimental re- 
sults as well as other lattice calculations, employ- 
ing a lattice discretisation of the Dirac action for 
the charm quark. To obtain reasonable results 
for charm quarks is crucial with respect to the 
calculation of the form factors for semi-leptonic 
B -^ D decays. Recent results for the semi- 
leptonic form factors have been presented as well. 
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